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We consider the tails of probability density function (PDF) for the velocity that satisfies Burgers 
equation driven by a Gaussian large-scale force. The saddle-point approximation is employed in the 
path integral so that the calculation of the PDF tails boils down to finding the special field-force 
configuration (instanton) that realizes the extremum of probability. For the PDFs of velocity and 



it's derivatives m'*' 



the general formula is found: \nV{\u^''''\) oc 



,(fc) 



|/Re'=) 



3/(fc+l) 



The forced Burgers equation 

dtu + udxU — vd^M = (/) 



(1) 



describes the evolution of weak one-dimensional acous- 
tic perturbations in the reference frame moving with the 
sound velocity [0. It is natural to assume the exter- 
nal force (\) to be (S-correlated in time in that frame: 
((^(ii,a;i)(/)(t2,a;2)) = i{t\ - t2)xixi - X2). Then the 
statistics of (j) is Gaussian and is completely character- 
ized by X. We are interested in turbulence excited by 
a large-scale pumping with some correlation length L 
so that X does not essentially change at a; < L and 
goes to zero where x > L. Besides L, the correlation 
function x niay be characterized by the parameter uj = 
[— (l/2)x"(0)]^/"^ having the dimensionality of frequency. 
Then e.g. x(0) ^ L?uj^. Developed turbulence corre- 
sponds to a large Reynolds number Re = Li^uj/v 3> 1. 

The physical picture of Burgers turbulence is quite 
clear: arbitrary localized perturbation evolves into shock 
wave with the viscous width of the front, that gives fc~'^ 
for the energy spectrum at Re ^ kL ^ 1 The 
presence of shocks leads to a strong intermittency, PDF 
of velocity gradients is substantially non-Gaussian |^ and 
there is an extreme anomalous scaling for the moments of 
velocity differences w = u{p) — u{—p): {w") oc (p/L) for 
n > 1 1^. Simplicity of the equation and transparency 
of underlying physics make it reasonable to hope that 
consistent formalism for the description of intermittency 
could be developed starting from Burgers equation . 
The most striking manifestation of intermittency is the 
deviation of the PDF tails from Gaussian. It has been 
realized recently that those tails can be found by consid- 
ering the saddle-point configurations (we call them in- 
stantons) in the path integral determining the PDF 
Note the simultaneous revival of the method of optimal 
fluctuation in the description of rare events in disordered 
metals P,p^. The instanton formalism has been applied 
to the Burgers equation first by Gurarie and Migdal pd| ] 
who found the right tail ln7'(u') ~ ~[w/{lljp)]'^ deter- 
mined by inviscid behavior of smooth ramps between the 
shocks. That cubic right tail has been earlier predicted 



by Polyakov from the conjecture on the operator prod- 
uct expansion it corresponds to the same right tail 
for the gradients lnV{u') ~ —{u'/loY derived by Go- 
toh and Kraichnan The PDF's Viu') and Viw) are 
not symmetric, the asymmetry is due to the simple fact 
that positive gradients are smeared while the steepening 
of negative gradients could be stopped only by viscosity. 
Here we describe the viscous instantons that give the left 
tails of PDF's and single-point velocity PDF. 

Even though some calculations are lengthy, the simple 
picture appears as a result. Since white forcing pumps 
velocity by the law vo^ oc while the typical time of 
growth is restricted by the breaking time t ^ L/w then 
the Gaussianity of the forcing lnV{(j>) oc — 0^/x(O) leads 
to InViw) ~ -[\w\/{Luj)]^. At a shock, ~ -vu' 
so that l-aV{u') ~ -[-u' / (ujRe)]^^'^ . These simple esti- 
mates are confirmed below by consistent calculations. 

Following we write the high-order moments of the 
velocity derivatives u^''^ — d^u (including velocity for 
/c = 0) and of the difference w as the path integrals: 

([uC^)]") =y"puPpexp|iZ + nln[u('=)(0,0)]}, (2) 

(w") = / PuPpexpfiZ-f nln[u(0,p) -u(0,-p)]}. (3) 



Here, I ~ j dt C{u,p) is the effective action with the 
Lagrangian C determined by the equation (n^) M-nffl: 



C — J dx {pdtU + pudxU - vpd^u) ^ l^j dxidx2PiXi2P2 ■ 

The main idea implemented here is that the high-order 
moments (for n 3> 1) are determined by the saddle-point 
configurations of the path integrals (||,||). The corre- 
sponding saddle-point equations are 



dtu + udxU — vd^u = —i / dx' x{x ~ x')p{x')^ (4) 



dtP + udxP + vdlp = 5{t)\{x). 



(5) 



We should substitute here A — [in/a\5'{x) with a{t) ~ 
dxu{t, 0) for gradients, A — —{in/w)[5{x — p) — 5{x + p)] 
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for differences etc. A solution of determines the 

moments in the steepest descent approximation: 



([uC-^)]") - [uC^) (0, 0)]"e'^--, (w-^) - [u;(0)]"e 



(6) 



An instanton solution describes the field-force configura- 
tion that corresponds to the rare fluctuation giving the 
main contribution into the high-order moment. This con- 
figuration can be also called optimal fluctuation. 

The initial conditions for the instanton equations were 
formulated in it ^ at < ^ — oo and p = at i > 0. 
The role of the last term in (||) is then reduced to the 
final condition p{x) = — A(a;) imposed on p at t = —0. 
Viscosity smears p if to move backwards in time so that 
at large negative time both fields u and p are zero. We 
will refer to that stage as vacuum since it is realized at 
the absence of the source A in the right-hand side of (||). 

In this letter, we describe general properties of instan- 
tons and find their dependence on n, which gives the 
form of PDF tails. Complete analytic descriptions of the 
instanton solutions will be published elsewhere. 

Since the Lagrangian does not explicitly depend on 
time, then the "energy" E is conserved by (0||): 



E = i dx{pudxU + vdxpdxu) 



1 



dxidx2PlXl2P2 



Since C does not explicitly depend on coordinates then 
the "momentum" J is conserved as well: iJ — J dxpdxU. 
Because of the conservation laws we should treat solu- 
tions of (^,||) with E = J = since they are zero in the 
vacuum. That gives the following saddle-point value of 
the effective action in (^: 2cxtr = / dtdxpdtu. 

Conservation laws help to understand general prop- 
erties of the solutions. We consider i = 0, substitute 
p{x) = — A(a;) into E and analyze the balance of different 
terms. For the gradient, E — — na(0) + oj^n^ / a?'{Q) — 
nvdxU{Q,Q) / a{Q) = 0. Difference between the cases 
of positive and negative a is now clearly seen. For 
a(0) > 0, the viscous contribution to the energy is 
unessential and two first terms can compensate each 
other (see below). Contrary, the instanton that gives 
a(0) < cannot exist without viscosity. For the ve- 
locity, J = gives dxu{Q, 0) = which leads to £' = 
x(0)n^/u^(0) - vdlu{{),Q) = 0. Without viscous term, 
energy conservation cannot be satisfied. Let us empha- 
size that the answer we shall obtain for the velocity PDF 
does not contain viscosity while it's consistent derivation 
requires the account of the viscous terms in the equations. 

Let us first describe the essentially inviscid instantons 
producing the right tails of the PDFs for gradients and 
differences |[|llHl|l. At t = 0, the field p is localized near 
the origin. At moving backwards in time the viscosity 
will spread the field p. Nevertheless, a positive velocity 
slope "compresses" the field p so that one can expect that 
the width of p remains much smaller than L. Then, it is 
possible to formulate the closed system of equations for 



the quantities a{t) and c{t) = —i J dxxp{t,x) since for 
narrow p and small x we can put J dx'x{x — x')p{t, x') — > 
—idxx{x)c{t) K, 2i(jj^xc{t): 



dtc = 2ac, dta 



(7) 



The instanton is a separatrix solution of (^. The initial 
condition a(0)c(0) = n by virtue of the energy conser- 
vation gives a(0) = uj^c^{Q)/n = wn^/'^. For differences, 
w = 2a(0)p. One can check that loxtr = * Jdtcdta ~ 
a(0)c(0) = n which is negligible in comparison with 
nln[a(0)] so that ((w')") ~ HO)]" ^ Lu^'n"'^^ which gives 
the right cubic tails of the PDFs InViu') ~ -{u'/uj)^ |l|] 
and \nV{w) ~ —[w/{pLo)]'^ Ml^- One can show that the 
width of p is much less than L through the time of evolu- 
tion T ^ n^^^^ijj~^ giving the main contribution into the 
action |l^. The right tails of V{u') and Viw) are thus 
universal i.e. independent of the large-scale properties of 
the pumping. Above consideration does not imply that 
the instanton is completely inviscid, it may well have vis- 
cous shock at a; ~ L, this has no influence on the instan- 
ton answer (since p is narrow) while may influence the 
fluctuation contribution i.e. predexponent in the PDF. 

The main subject of this paper is the analysis of the 
instantons that give the tails of V{u) and the left tails 
oiV{u') and V{'w) corresponding to negative a, w. Even 
though the field p is narrow at t — 0, we cannot use the 
simple system (^ to describe those instantons. The rea- 
son is that sweeping by a negative velocity slope provides 
for stretching (rather than compression) of the field p at 
moving backwards in time. As a result, the support of 
p{x) stretches up to L so that one has to account for 
the given form of the pumping correlation function xi^) 
a,t x ^ L. This leads to a nonuniversality of 'P{u) and 
of the left tails of V{u') and V{'w) which depend on the 
large-scale properties of the pumping. As we shall see, 
the form of the tails is universal, nonuniversality is re- 
lated to a single constant in PDF. Additional complica- 
tion in analytical description is due to the shock forming 
from negative slope near the origin. The shock cannot be 
described in terms of the inviscid equations so that we 
should use the complete system (^j^) to describe what 
can be called viscous instantons. 

Apart from a narrow front near x — Q, the velocity field 
has L as the only characteristic scale of change. The life 
time T of the instanton is then determined by the mo- 
ment when the position of p maximum reaches L due to 
sweeping by the velocity wo: T ^ L/uq. Such a veloc- 
ity uq itself has been created during the time T by the 
forcing so that uq ~ \c\maxTLu!^. To estimate the maxi- 
mal value of |c(i)|, let us consider the backward evolution 
from t — 0. We first notice that the width of p (which 
was zero at t = 0) is getting larger than the width of 
the velocity front ~ uo/a already after the short time 
~ a~^. After that time, the values of c and a are of 
order of their values at t = 0. Then, one may consider 
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that p{t, x) propagates (backwards in time) in the almost 
homogeneous velocity field uq so that 



V{u') oc exp[-C(-u7a;Re)3/2] 



(11) 



dxxupx ~ 2mo / dxp . 



dfC = -? 



The (approximate) integral of motion i J dxp can be es- 
timated by it's value at t = which is n/2uo. There- 
fore, we get Cmax — nT so that T ~ n~^^^u)~^ and 
Mo — Lujn^^^. At the viscosity-balanced shock, the ve- 
locity Uq and the gradient a are related by Ug ~ i/a so 
that a(0) ~ wRen2/3. 

Let us briefly describe now the consistent analytic pro- 
cedure of the derivation of the function c{t) that con- 
firms above estimates. We use the Cole-Hopf substitu- 
tion ^ for the velocity = —u'^/2v and introduce 
P — ivdxp/"^ ■ The saddle-point equations for \1/ and P 

a** - vdl^l + vF-^ ^ 0, (8) 
dtP + vdlP - vFP - 2v\' {x)5{t)^-^ = (9) 

contain F determined by dxF{t,x) = —i Jdx'xi^ ~ 
x')p{t,x')/2v'^ and fixed by the condition F(i, 0) — 0. 
We introduce the evolution operator lJ{t) which satisfies 
the equation dtU = HU with H{t) = v{dl - F). It is 
remarkable that one can develop the closed description 
in terms of two operators A = U~^xU and B = U~^dxU- 

dtA = -2iyB, dtl3 ^ -vFx{t,A) . 

Since we are interesting in the time interval when p{t, x) 
is narrow, it is enough for our purpose to consider x L 
where F{t, x) — c{t)x'^ijj'^ /2v'^ . Further simplification can 
be achieved in this case and the closed ODE for c{t) can 
be derived after some manipulations: 

{dtcf = Au^c^ + 16^2 + 4w^f? , 

where ^i — i J dxX{x)x and 4^2 = ^* JdxX{x)dx [xu{0, x)] . 
Integrating we get 



1 

t = - 



dx 



(0) Vw3x3 + 4^1 + 



(10) 



which describes c(t) in an implicit form. Further anal- 
ysis depends on the case considered. For the gradients, 
we substitute = n/ao and ^2 — -~n/2 and see that, 
as time goes backwards, negative c{t) initially decreases 
by the law c{t) = c(0) + 2nt until T = uj-'^{n/2)-^/^ 
then it grows and the approximation looses validity when 
c{t) approaches zero and the account of the pumping 
form x{^) at a; ~ L is necessary. Requiring the width 
of p{x) at this time to be of order L we get the esti- 
mate a(0) ~ ijjReri^/3 and thus confirm the above pic- 
ture. The main contribution to the saddle-point value 
(^ is again provided by the term [dxu{0,0)Y^ and we 
find {{u')") ~ [a(0)]" ~ (wRe)"^^"/^^ which corresponds 
to the following left tail of PDF at u' 3> a; Re 



For higher derivatives u^'^', by using ([lo| ) we get initial 
growth c{t) = c(0) + n{k + l)t which gives u('')(0,0) ~ 
^fc+i^i_fc^j^gfc leading to ([u^^)]") - ojRe^ L^-^n^^+^^l^ 
which can be rewritten in terms of PDF: 



V 



ocexp 



3/(fc-l-l) 



(12) 



Note that the non-Gaussianity increases with increasing 
k. On the other hand, the higher k the more distant is 



the validity region of (O): uC^) > u 



,(k) 



For the differences, ^1 ~ 2npQ/w and 4^2 = —n[l + 
2/9oUi(0, Pa)/w\ and we get (w") ~ {LwY^n'^^'^ which cor- 
responds to the cubic left tail 



V{w) oc exp{-B[w/(Ltj)]3} 



(13) 



valid at w 3> Luj. In the intermediate region Luj ^ w ^ 
puj, there should be a power asymptotics which is the 
subject of current debate |p|, p^JT7| . It is natural that p- 
dependence of V{'w) cannot be found in a saddle-point 
approximation; as a predexponent, it can be obtained 
only at the next step by calculating the contribution of 
fluctuations around the instanton solution. This is con- 
sistent with the known fact that the scaling exponent is 
ri- independent for n > 1: {w^{p)) oc p. 

For the velocity, A(a;) = — in(5(x)/u(0, 0) is an even 
function so that is a linear (rather than quadratic) 
function of x for narrow p: F{x) = x(0)&x/2i/^ with 
b — —i J dxp{x). Direct calculation shows that energy 
and momentum conservation makes h time independent: 
b — n/u{Q, 0). It is easy then to get the n-dependence of 
u(0,0): Velocity stretches the field p so that the width 
of p reaches L at T •2^ L/u{Q,0) while the velocity it- 
self is produced by the pumping during the same time: 
m(0,0) ~ x(0)6r = x(0)nr/2M(0,0) ~ nx(0)L/u(0,0). 
That gives ^(0, 0) ~ Luj-n}/^ and 

V{u) o^eyiY>{-D[u/{Luj)f} . 

The product Luj plays the role of the root-mean square 
velocity Ur,ns- The numerical factors C, B and D are de- 
termined by the evolution at i ~ T i.e. by the behavior 
of pumping correlation fmiction x{^) at x ~ L. 

We thus found the main exponential factors in the PDF 
tails. Complete description of the tails requires the anal- 
ysis of the fluctuations around the instanton which will 
be the subject of future detailed publications. Here, we 
briefly outline some important steps of this analysis. The 
account of the fluctuations in the Gaussian approxima- 
tion is straightforward and leads to the shift of Xcxtr in- 
significant at 71 ^ 1. However, the terms of the perturba- 
tion theory with respect to the interaction of fluctuations 
are infrared divergent (proportional to the total observa- 
tion time). That means that there is a soft mode which is 
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to be taken into account exactly. Such an approach has 
been already developed in [|j for the simpler problem of 
the PDF tails for a passive scalar advected by a large- 
scale velocity where the comparison with the exactly solv- 
able limits was possible. A soft mode usually corresponds 
to a global symmetry with a continuous group: if one 
allows the slow spatio-temporal variations of the param- 
eters of the transformation then small variations of the 
action appears. Our instantons break Galilean invariance 
so that the respective Goldstone mode has to be taken 
into account. Namely, under the transformation 



X — r, u{x) — > u{x — r) + w, r 



v{T)dT , (14) 



the action is transformed as X ^ X — i J dxdtpdtv. The 
source term J dxdtXu is invariant with respect to (|lj) 
for antisymmetric A(a;). To integrate exactly along the 
direction specified by ( p| ) in the functional space we use 
Faddeev-Popov trick inserting the additional factor 



1 = J Vv{t)S u(^t, 



v{T)dT ) — v{t) 



J . (15) 



into the integrand in (0,3[). Jacobian J is determined 
by a regularization of (04 ) according to our choice of the 
retarded regularization for the initial integral: at dis- 
cretizing time we put dtu + udxU — > (u„ — M„_i)/e -t- 
M„_iu^_2 (otherwise, some additional u-dependent term 
appears Q). The discrete version of ([l4[), p„(a;) 

un{x) 



Unix - e^. Vj) 



UNix) — vn gives 



dtu' it, / v{T)dT 



J = exp 



Substituting ( |l5| ) into (^,||) and making (|lj) we calcu- 
late / Vv as a Fourier integral (the saddle-point method 
is evidently inapplicable to such an integration) and con- 
clude that after the integration over the mode (|l^) the 
measure 'Du'D'pe^^ acquires the additional factor 



d'l'pit, x)dx 



6[u{t, 0)] exp 



u'{t,Q)dt 



-T 



The last (jacobian) term here exactly corresponds to the 
term u'Viu') in the equation for V{u') derived in |p|,p^. 
This term makes the perturbation theory for the fluc- 
tuations around the instanton to be free from infrared 
divergences, the details will be published elsewhere. 

Let us summarize. At smooth almost inviscid ramps, 
velocity differences and gradients are positive and lin- 
early related w(p) ~ 2pu' so that the right tails of PDFs 
have the same cubic form [^I JlllJl^ . Those tails are uni- 
versal i.e. they are determined by a single characteris- 
tics of the pumping correlation function x(r), namely, 



by it's second derivative at zero w = [— (l/2)x"(0)]^/'^. 
Contrary, the left tails found here contain nonuniversal 
constant which depends on a large-scale behavior of the 
pumping. The left tails come from shock fronts where 
ui^ ~ —lyu' so that cubic tail for velocity differences (^ 



corresponds to semi-cubic tail for gradients (|ll|). The 
formula (13) is valid for w ^ u^ms — Llo where 'P(u') 
should coincide with a single-point V{u) since the prob- 
ability is small for both u{p) and u{—p) being large si- 
multaneously. Indeed, we saw that the tails of \nV{u) 
at M ^ Urms are cubic as well. Note that ([l^) is the 
same as obtained for decaying turbulence with white (in 
space) initial conditions by a similar method employing 
the saddle-point approximation in the path integral with 
time as large parameter [Q. That, probably, means that 
white-in-time forcing corresponds to white-in-space ini- 
tial conditions. Note that if the pumping has a finite 
correlation time r then our results, strictly speaking, are 
valid for u,w L/t and u' <C l/r. 
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